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ABSTRACT 


This thesis examines the performance of single-carrier systems in a fading 
channel environment. A transmission through a fading channel results in changes in 
phase and amplitude that degrade system performance. A pilot symbol in a single-carrier 
system can determine the channel characteristics, allowing reliable system performance. 
The differential decision feedback and double differential-decision feedback algorithms 
compensate for Doppler error by updating the channel phase estimation between each 
symbol. 

To determine the performance of single carrier systems in a fading environment, 
this thesis used MATLAB simulations. Each algorithm was tested for a variety of 
conditions to determine performance characteristics. Simulations for QPSK, 16-QAM 
and 64-QAM were implemented to test bandwidth efficient modulation using pilot 
symbols. Higher order QAM signals are less tolerant to phase errors since the signal 
constellation points are more densely located. These simulations showed that both 
algorithms can effectively compensate for Doppler errors in a fading channel with similar 
performance. In most situations, the differential-decision feedback algorithm performed 


marginally better than the double-differential decision feedback algorithm. 
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EXECUTIVE SUMMARY 


As the use of digital wireless communications systems grows, more efficient 
modulation techniques are desired. Fading or multipath channels present additional 
challenges to designing wireless systems. Fading channel Doppler effects can induce 


phase errors, complicating receiver design. 


This thesis examines the feasibility of using pilot symbols to correct for phase 
errors in fading channels. Two proposed methods for recovering the phase information 
are examined. Both methods use a pilot symbol at fixed intervals to initially determine 
the channel’s phase error. Differential-decision feedback estimates the current channel 
phase shift from the last received symbol. This estimation is then used to decode the next 
symbol. Double-differential decision feedback extrapolates the last two received 


symbols to estimate the phase for the next symbol. 


To determine the effectiveness of the proposed Doppler correction techniques, 
MATLAB was used to model the techniques in a fading channel environment. M-ary 
phase-shift keying (MSFK) and M-ary quadrature amplitude modulation (MQAM) 
simulations were developed to test Doppler correction for bandwidth efficient modulation 
techniques. Transmitted signal power, measured in terms of E,/No, and the pilot symbol 
interval were varied to determine their effect on bit error ratios. Results from these 


simulations are compared to theoretical models for fading channels. 


The results from the simulations show that using a pilot symbol and then 
progressively estimating the channel phase error is an acceptable method. Bit error ratios 
on the order of 10° were achievable with the pilot symbols spaced 50 symbols apart. The 
simulations showed that the performance of the system improved as the pilot symbol 
interval decreased. Using the double-differential decision feedback algorithm added no 
noticeable performance improvement over the differential-decision feedback algorithm 


despite its use of two symbols to predict the future phase shift. 
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I. INTRODUCTION 


A. OVERVIEW 


Digital communications play an ever-expanding role in all aspects of modern life. 
As the information carried over wireless networks increases, better ways are needed to 
more efficiently make use of the bandwidth available. Complicating matters further, 
most real-world channels do not follow the basic Gaussian noise model. Instead, wireless 


transmissions arrive at a receiver from multiple paths. 


Communications transmissions over multiple paths are also known as fading 
channel communications. Often, the channel is constantly changing due to changes in the 
environment. If either the transmitter or the receiver is moving relative to the other, then 
it further complicates the problem. The multipath nature of the fading channel usually 
results in a received signal with varying phase and power. To compensate for the 
unknown conditions, the communications system must determine current channel 


characteristics. 


B. OBJECTIVE 


The primary goal of this thesis was to investigate the feasibility of using pilot 
symbols with single-carrier transmissions in a fading channel environment. A second 
goal was to investigate performance of single carrier transmissions in a fading channel 
when received signals were Doppler corrected using pilot symbol phase estimates. A 
third goal was to compare the performance of two pilot phase estimators. The MATLAB 
software program is used to simulate transmissions over fading channels. To reduce the 
frequency of pilot symbols, the estimated channel phase shift is updated as each symbol 
is received and decoded. This represents a new way to overcome Doppler shift induced 
phase errors in fading channels and, to the author’s knowledge has not been previously 


investigated. The simulations demonstrate the performance obtained when using pilot 


symbols for a variety of conditions using quadrature phase-shift keying (QPSK), 16- 
quadrature amplitude modulation (QAM) and 64-QAM. These represent a range of 
higher-order modulation techniques to demonstrate performance differences as 


bandwidth efficiency increases. 


Simulation of the signals for a variety of different configurations and channel 
characteristics shows how performance relates to various pilot symbol interval rates. 
Phase estimation using either one or two symbols to determine phase changes are 


compared to determine if there is a performance gain. 


c. OUTLINE 


This thesis consists of four chapters. Chapter I is a basic overview of goals and 
organization. Chapter II provides a theoretical background to single-carrier systems in a 
fading channel. Chapter III covers the MATLAB simulation implementation and results. 


Chapter IV reviews the overall results and their implications. 


UH. DOPPLER CORRECTION ALGORITHMS-A REVIEW 


A. FADING CHANNEL OVERVIEW 


Multipath channels present various challenges to wireless communication 
systems. In line-of-site (LOS) communications, the signal arrives at the receiver via only 
the direct-most path. A multipath signal arrives at the receiver from multiple angles and 
at slightly different times. Since signals take different paths to reach the receiver, they 


potentially interfere with each other creating changes in signal strength and phase. 


Fading channels are typically characterized by two parameters: coherence time 7, 
and coherence bandwidth B,. The coherence time is inversely proportional to the 
maximum Doppler spread, which is related to the channel’s varying phase and amplitude 
characteristics. When the coherence time is much larger than the symbol interval, the 
channel is described as slowly fading. This implies that the channel is static over at least 


one symbol transmission [1]. 


The coherence bandwidth is approximately the inverse of the multipath delay 
spread, the time difference between the first and last received path receptions. If the 
coherence bandwidth exceeds the signal bandwidth, the channel is a flat fading or 
frequency-nonselective channel. Channels where the coherence bandwidth is less than 


the signal bandwidth are frequency-selective channels [3]. 


There are many statistical models for fading channels, including Rayleigh, Rician, 
and Nakagami-m. Rayleigh fading models communication channels where there are 
many scattering objects to reflect signals with no LOS path. The amplitude for the signal 
after passing through a Rayleigh channel is modeled by the Rayleigh distribution 


f (x) =2xe™ (1) 
while the phase is uniformly distributed on the interval [- 2, 2] [1]. 
B. DOPPLER ERROR PERFORMANCE DEGRADATION 


With higher order modulation techniques, increased receiver phase errors lead to 


increased bit error rates since they make use of the phase information to differentiate 


3 


symbols. As the order of coherent modulation techniques increases, the effect of even 
small phase errors is magnified. Since 64-QAM symbols have less phase difference 
between them than 16-QAM or QPSK, it follows that for the same phase error, 64- 


QAM’s performance degrades more. 


For M-ary I-Q modulation schemes, for performance is given by [1] 


N Os 
P, ~ n QO min (2) 
log, M 20 


where dmin is the minimum Euclidian distance between two constellation points, NV, is the 





average number of nearest neighbors, and o” is the noise variance. When a phase error is 
introduced at the receiver, the distance between the received symbol and the expected 
decision vector decreases. The effective minimum distance dmin(€) as a function of phase 


error € is given by [1] 











di, (€) = min [Ts.-s,| “Its. -sf (3) 
min ij s, -s,| 
where 
T- cosé —siné Fi 
(sine | ) 


and s; is the I-Q vector location of the i" constellation point for all values of i and j that 


span the entire constellation set, i 4/. 


Using standard rectangular constellations where dmin is normalized to the square 
root of the energy per symbol (,/E, ) and with MATLAB to numerically solve Equation 


(3), we see in Figure | how quickly phase errors degrade performance for QPSK, 16- 
QAM, and 64-QAM. When the minimum distance reaches zero, the receiver is more 
likely to decode the received symbol incorrectly. This represents the limit in phase error 
that the receiver can tolerate. In the case of 64-QAM, any phase error exceeding 7.8° 
causes the receiver to become unreliable. Comparatively, QPSK performs much better 
when faced with phase error. With QPSK, at 7.8° dinin, has only decreased by a factor of 
0.85, which is only 1.36 dB in performance degradation. 
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Figure 1. Minimum normalized Euclidean distances as a function of phase error (dmin 


normalized to ,/E, ). 


Further examination of the results for QPSK demonstrates how significant the 
decreasing Euclidean distance is to performance with relatively small phase errors. 
Figure 2 shows QPSK performance with different induced phase errors in an additive 
white Gaussian noise (AWGN) channel. For a given energy per bit-to-noise power 


spectral density ratio E,/N,, theoretical performance was calculated using the special 


case for QPSK [1]: 





1 E, _ tf ; 
nio=to| a (cosé ine) +l 2— (osesine) (5) 


0 0 
As the phase error increases to 20°, the fall-off in performance quickly becomes 


significant. The performance decrease in Figure 2 demonstrates how the more-tolerant 


QPSK is susceptible to even small phase errors. Higher order modulation techniques, 
including 16-QAM and 64-QAM, experience similar performance degradations at smaller 


angles. 
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Figure 2. Effect of phase error on QPSK performance in AWGN. 


Small phase errors significantly degrade the performance of higher order 
modulation techniques. To maximize data throughput and increase reliability, the 


receiver must avoid large phase errors by updating its estimated phase regularly. 


PILOT SYMBOL DECISION FEEDBACK DEMODULATION: 
BACKGROUND 


In a fading channel, the received signal is subject to constant phase changes. A 
pilot symbol transmitted with the signal can be used to determine the current phase. 
Since the pilot symbol’s phase is known to both the transmitter and the receiver, the 
receiver can measure the difference to estimate the channel’s phase shift. In a single- 


carrier system, the pilot symbol should be sent periodically spaced between data symbols. 


Periodic updates determine the channel phase error at the time during the 
transmission of the pilot symbol; however, the channel phase can drift between pilot 
updates. _ Uncompensated phase changes between pilot symbols further reduce 
performance. Figure 3 shows basic operation of a single-carrier system using pilot 
symbols. A predetermined pilot symbol is multiplexed into the data at the transmitter at a 
fixed rate. The receiver separates the pilot symbol and measures the phase difference 
between the known reference and the received symbol to determine the channel’s phase 


error. This phase error is fed back to the receiver, which demodulates the signal. 


Complex Mapper Symbol- 
Matched Demulti- and to-bit 
Filter plexer Dectector converter 





Figure 3. Complex pilot symbol demodulation block diagram. 


Even in slowly fading channels, phase errors can quickly accumulate between 
pilot symbols. The change in phase at the receiver for a given Doppler shift f, and 
symbol period 7, after / symbols is given by [1] 

W = 2a fT (6) 


which shows that the phase error increases linearly between pilot symbol transmissions. 


1. Differential-decision Feedback Algorithm [1] 


The differential-decision feedback (D-DF) algorithm minimizes bit error rates by 
updating estimated phase error between pilot symbols. To compensate for the phase drift 
between pilot symbols, the estimated phase error can be updated continuously using data 
symbols. Using the pilot symbol, the receiver first calculates the current channel phase 
error. This phase correction is applied to the reception of the next symbol, reducing the 


phase error for that reception. 


To estimate the phase error after reception, the received value is compared with 
the decoded symbol value. Calculating the phase difference between the decoded and 
received value, we obtain an estimate of the current channel phase. The reception of the 
next symbol only requires the channel phase estimation from the previous symbol. The 
ha symbol uses the estimation from the (1-1)" for phase estimation. The phase error 


estimation can be calculated using [1] 
ee 3, +é, (7) 
where a, is the estimated pilot phase, and é, ,is the estimated phase difference for the 


previous reception This repeats until the next pilot symbol is transmitted. Figure 4 


shows a block diagram with the implementation of differential-decision feedback. 


Complex Mapper Symbol-— 
Matched Demulti- and to-bit 
Filter plexer Dectector converter 
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Figure 4. Complex pilot symbol demodulation with decision feedback block diagram. 


Correctly estimating the phase requires that no error was made when 
demodulating the symbol. If the symbol is decoded incorrectly, then the phase error 
calculation will be off. This erroneous phase can then lead to further decoding errors. In 
an attempt to limit the effect of symbol errors, a limit is placed on the maximum inter- 
pilot symbol phase correction. Both the type of modulation and expected nature of the 


fading channel must be used to determine the optimal phase error threshold. 


Ze Double-differential Decision Feedback Algorithm [1] 


The double-differential decision feedback (DD-DF) algorithm is an attempt to 
improve the performance of the differential-decision feedback algorithm by using the 


phase error from the last two received symbols to predict what the phase error will be for 
8 


the next symbol. Since the channel’s phase is continually changing, the last symbol’s 
estimated phase lags the current phase and is a source of error. Use of the last two 
received symbols to estimate the symbol’s phase should improve the estimate, reducing 


the phase error and improving performance. 


Since the channel is assumed to be slow fading, the phase rate of change between 
two symbols should remain constant. This means that a simple linear regression can be 
used to accurately model phase changes. To estimate the phase on the next reception, the 
estimated phase change over the last two received symbols is extrapolated. The 


estimated phase for the pe reception can be calculated using [1] 
0 = 3, +2€,,— Ey (8) 
where a is the estimated pilot phase, and é,is the estimated phase difference for the i 


reception. 


D. CHAPTER SUMMARY 


This chapter reviewed the nature of fading channels and covered how the 
potential phase error caused by a fading channel degrades communication system 
performance. The differential and double-differential decision feedback algorithms were 
proposed to correct channel phase errors. The next chapter will show the results of 


simulations using these algorithms for fading channels. 
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10 


HiIl. PERFORMANCE OF MPSK AND MQAM IN A MOBILE 
FADING CHANNEL WITH DOPPLER CORRECTION 


A. SIMULATION OVERVIEW 


Computer simulation provides a useful tool to ascertain the performance of single 
carrier systems in a fading channel. Use of a model of the system and the channel 
permits the exploration of performance characteristics for a variety of conditions. The 
MATLAB software package was used to create a model that implemented both the 


differential and double-differential decisions feedback algorithms. 


To create the simulation, models for QPSK, 16-QAM, and 64-QAM were 
implemented. To verify that these models worked, they were initially tested with an 
AWGN channel to confirm that their performance matched well-known theoretical 
results. To reduce the number of computations required, the simulation modeled 
complex-valued baseband signals. Nyquist sampling requirements of a modulated signal 
in the passband greatly increases computational requirements and lengthens simulation 


run times. 


To model the modulation schemes, the modem object in the communications 
toolbox was used. Use of the modem class allowed the development of generic functions 
to implement both the D-DS and DD-DS algorithms. The modem class is capable of 
representing most types of PSK and QAM modulations. To reduce bit errors given a 
symbol error, grey coding was used. Once the modem modulator object has been 
created, the modulate function is used to encode a string of numbers. The input string 
is limited in values from 0 to M-1 and the output is a series of complex values. A modem 
demodulator object demodulated the baseband complex envelope signal using hard 


decisions. 


Simulation of the systems’ performance was accomplished using Monte Carlo 
trials. To determine bit error ratios, a test signal uniformly distributed on the discrete 
interval [0, M-1] was encoded and then passed through a simulated channel. The 


received channel was then demodulated and compared to transmitted signal to determine 
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the number of bit errors. To reduce simulation memory requirements, blocks of 10,000 
symbols were simulated at time. Since determining the bit error rate accurately requires 
many more samples, blocks were simulated repeatedly until enough trials had taken 
place. As each trial was run, a tally of error bits and total bits was maintained. The exact 
number of trials depended what the BER for that simulation was indicating with no fewer 


than 300 trials run. 

To determine the performance of the system in various conditions, simulations 
were performed for E, /N, values ranging from 10 to 61 dB at 3 dB intervals. Simulation 
over this range showed typical performance over a wide range of expected conditions. 


Results from the simulation were graphed for comparison and benefit analysis. The 


theoretical performance in Rayleigh fading with no phase error given for QPSK are [1] 


_! [_4/No 
a=z! oS] a 


E 
1.5(2i-1) log, M 
= JM /2 : 
naa ‘ 1 : No 
VM )\log,M } M -1+1.5(2i-1)’ > log, M 


0 





and for MQAM by [9] 











(10) 


To simulate the fading channel, the rayleighchan object from the 
communications toolbox was used. For frequency-nonselective channels, the properties 
of the channel are determined by two parameters, the sample rate and the Doppler 
frequency. Together these determine how fast the channel varies. The Jake’s Doppler 
power spectrum as a function of Doppler frequency shift v for a maximum Doppler shift 
fal3] 


0) ee 


d ioe ae 
a 


was used for the simulation since it has been determined to accurately model the Doppler 








VSS, (11) 


effect between mobile sources. To speed up simulation § results the 





ResetBeforeFiltering property was set so that each batch of 10,000 symbols 
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passed through an independent channel. To model the fading channel, the filter 
function was used with the rayleighchan object. After passing through the 


Rayleigh channel, Gaussian noise was added to the signal to model receiver noise levels. 


To implement the D-DF and DD-DF decoding algorithms, new functions were 
written. The highest energy symbol in the constellation was used for a pilot symbol and 
was inserted at given intervals in the input data before encoding. Use of the highest 
energy symbol for the pilot makes it less susceptible to deep fades and other sources of 
error. When decoding the received signal, the function steps through each symbol 
individually with the demodulate function to establish the received symbol. To 
determine the channel phase at the pilot symbol, the received complex symbol is divided 


by the expected pilot symbol, calculating the angular difference. 


Once the receiver has an estimate of the channel phase error from the pilot symbol, 
it can use that to correct Doppler errors in succeeding symbols. The next symbol is 
rotated to counter the estimated channel phase error. After decoding, the rotated symbol 
is compared to the decoded value’s constellation point. The angular difference is then 
used to update the next’s symbols estimated phase error. The process repeats until the 
next pilot symbol transmission. Since the simulation uses complex signaling, calculating 


the error is just complex division. 


Implementing the DD-DF algorithm uses the same basic steps as the D-DF. The 
first symbol after the pilot again is used to estimate the phase error for the second symbol. 
The difference in implementation comes after the reception of the second symbol after 
the pilot. Now the last two phase error estimates are used to predict the future phase 
errors. The last two corrections to the pilot phase error are retained to project the future 


Doppler error. 


Since there is a possibility that a symbol could be incorrectly demodulated, the 
algorithm must be able to reject spurious phase corrections. This is implemented with a 
threshold phase error that all estimates must fall below. If the threshold is exceeded it is 
assumed that the update is not constructive. The last good update is used to estimate the 


current phase error. Allowing bad phase estimates to be used greatly degrades 
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performance since an erroneous phase error leads to further demodulation errors which 


continue to propagate until the next pilot symbol transmission. 


Determination of the optimal phase error tolerances depends on the type of 
modulation and channel nature. The tolerance must be small enough to prevent decoding 
errors from creating spurious phase estimations errors. It also must be large enough so 
that it does not eliminate phase shifts that are a result of the Doppler shift. The upper 
limit on the tolerance depends on the modulation type since that determines the angular 
space between constellation points. The lower limit is determined by the channel 


Doppler shift and the signal symbol period. 


Using the results in Figure 1, we get a good basis to determine the upper tolerance 
limit. If the tolerance limit is chosen so that dmin is zero or less, then it is likely that a 
demodulation error due to noise will be interpreted as a large phase shift. A tolerance 
that corresponds to a din that is approximately 1.0 should provide ample protection from 
excessive phase changes. For QPSK, a tolerance of 0.35 radians was used in the 
simulations. In the 16-QAM simulations, the phase error correction was limited to 0.2 
radians. The 64-QAM tolerance was set at 0.1 radians. All of these tolerances proved 


sufficient. 


B. QPSK SIMULATION PERFORMANCE 


To test the pilot-aided single-carrier system, initially QPSK was simulated to 
verify that the algorithms were working as designed. A slowly fading channel where the 
system operated 1 Msps was initially used. The maximum Doppler shift was initially set 
at 50 Hz, a level that introduced small Doppler errors without pushing simulation limits. 
This model allowed verification that the phase error term correctly tracked the Doppler 


shift. Using Equation (6) to calculate the maximum phase shift per symbol, we get 


Y = 2nT, f,l = 2 (1/1108 / sec) (S0Hz) = 0.018" (12) 


This means that it would take at least 2500 symbols for the phase to drift to 45°, where 


dmin reaches zero. 
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To determine the optimum pilot interval, the first simulation held E,/N, constant 


at 40 dB. The bit error ratio (BER) was measured as the pilot interval was decreased 
from 5000 to 10 symbols. The results from this simulation are shown in Figure 5. 
Performance rapidly improved with decreasing pilot interval until the pilot interval 
reached a range from approximately 100 to 200 symbols. Significantly, this represents 
approximately 10% of the symbols required for the phase error to exceed the angle where 


dinin approaches zero. 
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Figure 5. QPSK probability of bit error vs. pilot symbol interval for E, /N, = 40 dB, R, 
= 1 Msps, and fa=50 Hz. 


Given that performance gains taper off after the pilot interval is reduced to 200, 
further simulations were run to investigate the algorithm’s performance. The next 
simulations examined how BER varied with different signal strengths with the pilot 
interval set at 50, 100 and 200. The results from these simulations are shown in Figure 6. 


All three cases closely followed the ideal case until E,/N, exceeds approximately 35 
15 


dB. Errors in the channel phase estimation cause the simulations to have poorer 


performance than the ideal. At E,/N, greater than approximately 40 dB, all three 


simulations show significantly poorer performance than the ideal. 
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Figure 6. QPSK D-DF probability of bit error vs. E,/N, for R; = 1 Msps and fa = 50 
Hz. 


After establishing the baseline performance of the D-DF algorithm in QPSK, the 
DD-DF algorithm was simulated to determine its effect on performance. The DD-DF 
performs similarly to the D-DF algorithm, as can be seen in Figure 7. There is no 
consistent improvement in performance. At lower E,/N, , using two symbols to 
estimate the future phase performed marginally worse. At higher E,/N,, the results for 


D-DF and DD-DF crossover, with D-DF becoming worse. 
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Figure 7. QPSK DD-DF probability of bit error vs. E,/N, for R;= 1 Msps, fa=50 Hz, 
and pilot interval = 100 symbols. 


Cc. 16-QAM SIMULATION PERFORMANCE 


Extending the simulations from QPSK to 16-QAM, we evaluated performance to 
determine the effectiveness of using a pilot symbol as well as using the D-DF and DD- 
DF algorithms. The simulation channel was again chosen to vary slowly over time using 
the same channel characteristics as before. The maximum Doppler phase change 
between symbols remained 0.018°. Even though the phase was changing at the same 
rate, it effects were more significant. For 16-QAM, dmin reaches zero at 17.3°, which 
occurs after at least 961 symbols. While significantly less than 2500, this is still 


significantly greater than the pilot symbol interval. 
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The results of the simulations for 16-QAM are shown in Figure 8. Both the 100 
and 200 symbol pilot interval simulations track the no-phase error performance to about 
40 dB. As expected, the 100 symbol interval simulation performs slightly better than the 


200 symbol interval simulation. 
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Figure 8. 16—-QAM probability of bit error vs. E,/N, for Rs = 1 Msps and f, = 50 Hz. 


D. 64-QAM SIMULATION PERFORMANCE 


The initial 64-QAM simulations closely follow the 16-QAM simulations. For the 
same channel, din reaches zero after 433 symbols if the phase is permitted to change 
with no updates. Here, the benefit to reducing pilot symbol interval is more pronounced 


for lower E,/N,. Reducing the pilot symbol interval from 100 symbols to 10 symbols, 


we see an improvement of 12 dB, with a corresponding decrease of throughput from 99% 
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to 90%. At E,/N, greater than 50 dB, the performance gains are even more pronounced. 


While the higher pilot interval results trend away from the ideal case, the 10 interval line 


continues to follow it. 
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Figure 9. 64-QAM probability of bit error vs. E,/N, for Rs = 1 Msps and fa= 50 Hz. 


E. PERFORMANCE IN ENVIRONMENTS WITH HIGHER DOPPLER 
SHIFTS 


While the previous results verify that both the D-DF and DD-DF algorithms 
effectively correct phase error, the environment simulated is one where the phase is not 
changing by a significant amount each symbol. This section examines the results for 
QPSK and 64-QAM in an environment with higher Doppler shifts. QPSK and 64-QAM 
were chosen since they represent extremes in phase error tolerance. For these simulations, 


the symbol rate was decreased and the Doppler shift was increased, causing a larger 
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maximum phase shift between each symbol. This section also more closely examines the 


relationship between the D-DF and DD-DF algorithms. 


1. QPSK 


To begin the study with higher Doppler rates, D-DF QPSK was simulated for 
channels with increasing maximum Doppler shift. The pilot interval was kept fixed at 
100 symbols to show the effect of the Doppler shift. The first channel was modeled with 
R, set at 200 ksps and fa set at 50 Hz. This increased the maximum phase shift between 
symbols calculated using Equation (6) to 0.09°. Even at five times the original phase 
change, this is still a small amount of phase change between pilot intervals. The 
performance in this case is not much worse than the previous case, as can be seen in 
Figure 10. Performance is closely matched until E,/N, is greater than approximately 30 
dB, where bit error ratios begin to greatly exceed the ideal case. 

Also shown in Figure 10, increasing the maximum Doppler frequency to 120 Hz 
and further decreasing the symbol rate to 100 ksps, we obtain a significant impact on 
performance. The potential phase error, from Equation (6), is 

W =2nT. fl =27(1/10° /sec)(120Hz) = 0.432° (13) 
In the time it takes to transmit 100 symbols, the phase error could potentially increase to 


43.2°. This much phase error would severely impact the system’s performance. 
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Figure 10. QPSK probability of bit error vs. E,/N, for pilot interval at 100 symbols. 


Further investigation of the channel with f7 =120 Hz and R, =100 ksps reveals that 
using post-estimated feedback allows for acceptable wireless bit error rates. Even with 
the interval maintained at 100 symbols, P;, is 0.002 at E,/N, = 45 dB. Without using 
phase estimates from received symbols to update the current phase, the performance 
cannot approach this level. The simulation results for QPSK using both the D-DS and 
DD-DS algorithms for pilot intervals at 100, 50, 25 and 10 symbol is shown in Figure 11. 
In all cases, the DD-DS algorithm performs just noticeably worse than D-DS algorithm. 
Performance for all six cases starts closely bunched and, as E,/N, improves, the 
simulation with closer spaced pilot symbols begins to perform better. However, even 
with only a 10 symbol interval, performance is significantly worse than the ideal case for 


E,/N, > 1. 
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Figure 11. QPSK probability of bit error vs. E,/N, for R; = 100 ksps 
and fy = 120 Hz. 


2. 64-QAM 


Using the same channel with R,; = 100 ksps and fa =120 Hz, we get even greater 
demands on phase error compensation for 64-QAM. Because of its tightly packed 
constellation 64-QAM, is much more susceptible to phase error than QPSK. After only 
7.8° of phase rotation due to error, performance becomes unreliable. Since with every 
symbol transmission the phase can change up to 0.432°, it takes only 18 symbols of 


uncorrected phase errors to exceed 7.8°. 


The results for the 64-QAM simulations through this channel are shown in Figure 
12. Both the D-DS and DD-DS algorithms were simulated with pilot symbols at intervals 
of 50, 25, 10 and 5. For all eight simulations, the performance rapidly approaches a floor 


and then remains almost constant for 64-QAM. This set of simulations indicates a few 
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differences between D-DS and DD-DS performance. At higher E,/N,, the DD-DS 
algorithm performs slightly better than D-DS algorithm. The difference is not much, but 
it is measurable, especially for the 50 symbol interval case. At low E,/N, , the 
simulations using pilot intervals of 10 and 25 symbols, the DD-DS algorithm performs 


considerably worse than the D-DS algorithm. 


Examination of both QPSK and 64-QAM in a situation where the phase is 
changing rapidly with respect to the pilot interval demonstrates the full capabilities and 


limitations of using pilot symbols with differential-phase correction in a fading channel. 
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Figure 12. 64-QAM probability of bit error vs. E,/N, for R; = 100 ksps 
and fy = 120 Hz. 
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F. CHAPTER SUMMARY 


This chapter examined the MATLAB simulation results for the differential and 
double-differential decision feedback algorithms. QPSK, 16-QAM, and 64-QAM 
simulation results were shown for a variety of fading environment. The next chapter will 


review the implications of these results. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 
1. Concept Validity 


Simulation of single-carrier systems where the signal is transmitted over a fading 
channel for a broad range of conditions has demonstrated the advantage of using pilot 
symbols for wireless communications in a mobile environment. Use of pilot symbols 
allows a narrowband single-carrier system to estimate the fading channel phase and then 
correct for it. This permits the use of higher order bandwidth-efficient modulation 


schemes. 


Both the differential and double-differential decision feedback algorithms 
estimate phase corrections, which permits the use of MQAM techniques where accurate 
phase reconstruction is required to demodulate the signal. The performance of QPSK in 
a fading channel where the coherence rate is as high as 0.1 times the pilot symbol 
transmission rate is shown in Table |. Under these conditions, the bit error ratios are 
greater than the ideal performance, but forward error correction can be used to bring the 
system performance down to acceptable levels. Reduction of the throughput to 75% 


greatly improves the bit error rates without significant loss in available data rate. 
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Simulation | Pilot symbol interval | Throughput | Data bit rate | BER 

D-DS 100 99% 198 kbps 0.00247 
D-DS 50 98% 196 kbps 0.00111 
D-DS 25 96% 192 kbps 0.00051 
D-DS 10 90% 180 kbps 0.00019 
D-DS a 80% 160 kbps 0.00008 
D-DS 4 75% 150 kbps 0.00006 
DD-DS 100 99% 198 kbps 0.00269 
DD-DS 50 98% 196 kbps 0.00125 
DD-DS 2) 96% 192 kbps 0.00060 
DD-DS 10 90% 180 kbps 0.00021 
DD-DS 5 80% 160 kbps 0.00008 
DD-DS 4 75% 150 kbps 0.00001 


Table 1. | Summary of QPSK bit error ratios compared to throughput 
with R; = 100 ksps, fa = 120 Hz, and E,/N, = 46 GB. 


The simulations results for 64-QAM using pilot symbols with phase error 
correction is summarized in Table 2, where the performance of 64-QAM for a slow 
fading channel when the coherence rate is significant compared pilot symbol rate is 
shown. Since 64-QAM is more susceptible to phase error than QPSK, the pilot symbol 
must be sent at a higher rate. The tradeoff for the high pilot symbol rate and lower 
throughput is better bit error ratios. With throughput reduced to 75%, the BER is further 
reduced, but on the order of 10” is still significant. Again, the channel bit error ratios are 


large enough to require forward error correction for reliable digital communications. 
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Simulation | Pilot symbol interval | Throughput | Data bit rate | BER 

D-DS 50 98% 588 kbps 0.04213 
D-DS ZS 96% 576 kbps 0.01228 
D-DS 10 90% 540 kbps 0.00249 
D-DS 2 80% 480 kbps 0.00077 
D-DS 4 75% 450 kbps 0.00053 
DD-DS 50 98% 588 kbps 0.03688 
DD-DS 25 96% 576 kbps 0.01252 
DD-DS 10 90% 540 kbps 0.00256 
DD-DS 5 80% 480 kbps 0.00077 
DD-DS 4 75% 450 kbps 0.00055 























Table 2. Summary of 64-QAM bit error ratios compared to throughput 
with R; = 100 ksps, fa = 120 Hz, and E,/N, = 46 GB. 


Despite the expectations for the double-differential decision feedback algorithm 
to improve phase error correction, it performed either similarly or worse to the 
differential decision feedback algorithm except for in a few rare cases. The only time 
where the DD-DF algorithm showed notable improvement was for high E,/N, , where 
performance appeared to reach an asymptotic limit. In most simulations, the performance 
of the DD-DF algorithm mirrored that of the D-DF algorithm. In some cases, noticeably 
at low E,/N,, the BER obtained using the DD-DF algorithm was significantly higher 


than that obtained for the D-DF algorithm. 

A possible explanation for the increased BER with the DD-DF algorithm is in 
how it extrapolates the future phase estimate from the last two phase estimates. Since the 
phase error estimates are noisy, extrapolation increases the effect of the noise. At low 


E,/N,. there is more noise, which explains why the trend is more noticeable there. A 


2d 


potential solution would be to use more than two symbols to predict the future phase. For 
short sequences, a linear regression could be used to estimate the next symbol’s phase 
error. If the rate of phase change is near constant over a few symbols periods, the use of 
more than just two symbols for linear extrapolation would provide a better solution than 


extrapolation from only two points. 


When evaluating the usefulness of using pilot symbols with phase error correction 
for slowly fading flat channels, both algorithms do impose some limitations on their use. 
The most noticeable limit comes from the tolerance threshold set to reject spurious phase 
estimates from incorrectly-demodulated symbols. While it may seem a good choice to 
increase the tolerance to handle faster fading channels, simulations with the tolerance set 
above the recommended values from Figure 1 showed considerably poorer performance. 
This effectively limits the maximum Doppler shift for a system based on the modulation 
type and symbol rate. Ifa larger tolerance is needed, a modulation technique with greater 


angular separation between the constellation points should be considered. 


Ze Summary 


The differential and double-differential decision feedback algorithms are a useful 
way to combat phase error due to Doppler shift in a fading channel using a single-carrier 
signal. The simulations indicate that acceptable channel bit error ratios are achievable. 
Instead of using bandwidth to continually transmit a pilot signal, pilot symbols can be 
sent at a rate appropriate for the channel. Since the pilot symbol is only sent with the 


data transmission, available bandwidth can be used to increase data transmission rates. 


B. RECOMMENDATIONS 


Further research into single carrier transmissions using pilot symbols in a fading 
environment should examine improved methods to estimate phase error. Using a history 
containing more than two symbols should allow more accurate predictions of phase error 
and improve performance. Investigations into the use of multiple symbols should 


include methods besides just linear extrapolation. 
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APPENDIX-MATLAB SOURCE CODE 


A. QAM16_d_min.m 





function x = QPSK_d_min(e) 
& Jeff Smith 
% NPS 
T = [cos(e), -sin(e); sin(e), cos(e)]; 
dmin = 1e9; 
for xh = —leZ2eL 
for yl = -1:2:1 
si = [xl; yll; 
for x2°=—1Le221 
for y2 =-1:2:1 
sj} = [x2; y2]; 
if (or (x1l~=x2,yl~=y2) ) 
d = (dist(sj',T*si) .*%2- 


dist (si',T*si).*2)/dist(si',sj); 
if (d < dmin) 
dmin = d; 
end 
end 
end 
end 
end 
end 
x = dmin; 


B. QAM16_d_min.m 


function x = QAM16_d_min(e) 
& Jeff Smith 


% NPS 
T = [cos(e), -sin(e); sin(e), cos(e)]; 
dmin = 1e9; 
for xl = =332<3 
for yl H=3:2223 
si = [xl; yl]; 
fOr XZ -H0- 32223 
for y2 =-3%233 
sj = [x2; y2]; 
if (or (x1l~=x2,yl~=y2) ) 
d = (dist(sj',T*si) .%2- 


dist (si',T*si).*2)/dist(si',sj); 
if (d < dmin) 
dmin = d; 
end 
end 
end 
end 
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end 
xX = dmin; 
C. QAM64_d_min.m 
function x = QAM64_d_min(e) 
SJeff Smith 
SNPS 
T = [cos(e), -sin(e); sin(e), cos(e)]; 
dmin = 1e9; 
for xl = -7:2:7 
for yl = -7:2:7 
si = [xl; yl]; 
for x2 = -7:2:7 
for y2 =-7:2:7 
sj = [x2; y2]; 
if (or (x1l~=x2,yl~=y2) ) 
d = (dist(sj',T*si). 
dist (si',T*si).*2)/dist(si',sj); 
if (d < dmin) 
dmin = d; 
end 
end 
end 
end 
end 
end 
x = Admin; 
D. Phase_e_plot.m 
SJeff Smith 
SNPS 
% plots dmin for QPSK, QAM16, QAM64 
for a = 1:100; 
P(a) = QPSK_d_min(a/100*pi/4) 
16(a) = QAM16_d_min(a/100*pi/4) 
64(a) = QAM64_d_min(a/100*pi/4) 
end 
figure; 
hold on; 
plot ([1:100]/100*45,QP, a 
plot ([{1:100] i Mowe ie '); 
plot ([1:100]/100*45,Q64, 'r'); 
axis([0,45,0,2]); 
legend ('QPSK', 'QAM16', 'QAM64"'); 
E. QPSK_p_e_t.m 


end 





sdJeff Smith 
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AQ- 

















SPlots theoritical phase error for QPSK 
c = 'bgrcmykb'; 
EDNo = 0:.5:16; 
hold on; 
for e = 0:4 
phi = e*5*pi/180; 
P_err = 1/4*erfc(sqrt (10.% (EbNo/10) ) .* (cos (phi) - 
sin(phi)))+1/4*erfc (sqrt (10.* (EbNo/10)) .* (cos (phi)+sin(phi))); 
semilogy (EbNo,P_err, [c(e+1l) '--']); 
hold on 
end 


F. D_DF_ demod.m 


function r = D_DF_demod(m, hmod, hdemod, pilotSymbol, 
pilotiInterval) 

& Jeff Smith 

%& NPS 

% Demodulates using D-DF algorithm 


LOLs ile 


ag 
t 


zeros (1,length(m)); 
zeros (1,length(m) ); 





zpilot = modulate (hmod, pilotSymbol)); 


for n = 1:length(m) 





if rem(n,pilotInterval) == 1 
r(n) = pilotSymbol; 
h =m(n) / zpilot; 
e = -l*angle(h); 
a = abs(h); 
f(n) = 6; 
else 
r(n) = demodulate (hdemod,m(n)/a*conj(h) /abs (h) ); 


z = modulate (hmod,r(n)); 
hl = m(n) / Zz; 


el = -angle(hl); 
if or(((abs(el-e))< Tol),and(abs(e)>(pi-Tol),abs (el) >(pi- 
Tol) )) 
e-=el; 
h = hil; 
end 


end 
end 


G. DD_DF_demod.m 


function r = DD_DF_demod(m, hmod, hdemod, pilotSymbol, 
pilotiInterval) 

SJeff Smith 

SNPS 

sUses DD-DF algorithm to demodulate 
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6 
| 


= zeros(1l,length(m)); 
= zeros(l1,length(m)); 





Fh 
| 


zpilot = modulate (hmod, pilotSymbol); 


for n = 1:length (m) 








if rem(n,pilotInterval) == 1 
r(n) = pilotSymbol; 
h = m(n) / zpilot; 
e = angle(h); 
el = 0; 
e2 = 0; 
a = abs(h); 
else 
r(n) = demodulate (hdemod,m(n) /a*conj(h) /abs (h) *exp (- 


z = modulate (hmod,r(n)); 








h = angle(m(n)/z); 
if or(((abs(eh-el))< Tol),and(abs (ehte) > (pi- 
Tol),abs(elte)>(pi-Tol))) 
e2 =el; 
el = eh; 
end 
end 
end 
D_DF_QPSK.m 


function Pb = D_DF_QPSK(EbNo, pilotInterval) 
% D_DF_QPFK 
& Jeff Smith 
%& NPS ECE 
%& April 23, 2008 

















Rs = le5; % symbol rate 
fl = 1*Rs; %& £1 

fs = f£1*1; % sample rate 
fd = 120; 

M = 4; % 


k = log(M)/log(2); 
testSize = 1e4; 





currentPerr = 0; 
err = []; 
NumTests = 0; 
totalErrors = 0; 


chan = rayleighchan(1/(fs),fd); 
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chan.StoreHistory = 0; 
chan.ResetBeforeFiltering = 1; 





hmod = 
modem.pskmod('M',M, 'PhaseOffset',pi/4,'SymbolOrder', 'Gray'); 
hdemod = modem.pskdemod (hmod) ; 





while and(or(totalErrors < 5000,NumTests < 400),NumTests <3500) 


NumTests = NumTests + 1; 
s = fix(rand(1,testSize) *M); % create random signal 
s(l:pilotInterval:testSize) = 3; 





m = modulate (hmod,s); 
m = rectpulse(m,fs/Rs) ; 
m = filter(chan,m)j; 


m = awgn(m,EbNo + 10*1logl10(Rs/fs)+ 10*1log10(k)); 








B- 


m = intdump(m, fs/Rs); 
r = D_DF_demod(m, hmod,hdemod,s(1),pilotInterval) ; 




















totalErrors = totalErrors + biterr(r,s); 
currentPerr = totalErrors/ (NumTests* (testSize- 
floor (testSize/pilotInterval))*k); 
end 
Pb = currentPerr; 


I. D_DF_QPSK_plot.m 


% Jeff Smith 
% NPS 
% Creates EbNo plots for QPSK 





eal 





bNo = 10:3:61; 





err = zeros(1,length(EbNo) ); 


for n = 1:length (EbNo) 














err(n) = D_DF_QPSK(EbNo (n),100) 
end 
semilogy (EbNo (1:length(err)),err) 
hold on; 
semilogy (EDNo, berfading(EbNo, 'psk',4,1),'--"'); 
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J. D_DF_QAM.m 























function Pb = D_DF_QAM(EbNo, pilotInterval) 

% D_DF_QAM 

& Jeff Smith 

%& NPS ECE 

%& April 23, 2008 

Rs = led; % symbol rate 

fl = 1*Rs; %& f1 

fs = f1*1; % sample rate 

fd = 120; 

M = 64; % 

k = log(M)/log(2); 

testSize = 10000; 

currentPerr = 0; 

NumTests = 0; 

totalErrors = 0; 

chan = rayleighchan(1/(fs),fd); 

chan.StoreHistory = 0; 

chan.ResetBeforeFiltering = 1; 

hmod = modem.qammod('M',M, 'SymbolOrder', 'Gray'); 

hdemod = modem.gqamdemod (hmod) ; 

z = 0:M-1; 

y = modulate (hmod, z); 

Pavg = sum(abs(y) .*2)/M; 

[a, b] = max(abs(y)); 

while and(or(totalErrors < 4000,NumTests < 300),NumTests <2500) 
NumTests = NumTests + 1; 
s = fix(rand(1,testSize) *M); % create random signal 
s(l:pilotInterval:testSize) = b-1; 
m = modulate (hmod,s)j; 
m = rectpulse(m,fs/Rs) ; 
m = filter(chan,m); 
m = awgn(m,EbNo + 10*1lo0g10(fs/Rs) +10*10g10 (k) 

,LO*logl10 (Pavg) ); 

m = intdump(m, fs/Rs); 
r = D_DF_demod(m, hmod, hdemod,s(1),pilotInterval) ; 
totalErrors = totalErrors + biterr(r,s); 
currentPerr = totalErrors/ (NumTests* (testSize- 


floor(test 
end 






































Size/pilotInterval))*k); 
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Pb = currentPerr; 


K. D_DF_QAM_plot.m 


SJeff Smith 
SNPS 
SPlots QAM BER curve 





EbNo = 10:3:61; 








err = zeros(1,length(EbNo) ); 





for n = 1:length(EbNo) 





err(n) = D_DF_QAM(EbNo(n),5) 


end 








semilogy (EbNo(1:length(err)),err) 
hold on; 
semilogy (EbNo, berfading(EbNo, 'qam',64,1),'--'); 














L. DD_DF_QPSK.m 


function Pb = DD_DF_OPSK(EbNo, pilotInterval) 
% DD_DF_QPSk..m 

& Jeff Smith 

S$ NPS ECE 
% April 23, 2008 

















Rs = le6; % symbol rate 
fl = 1*Rs; % £1 

fs = f1*1; % sample rate 
fd = 50; 

M = 4; % 

k = log(M)/log(2); 

testSize = 1le4; 

currentPerr = 0; 

NumTests = 0; 

totalErrors = 0; 








chan = rayleighchan(1/(fs),fd); 
chan.StoreHistory = 1; 
chan.ResetBeforeFiltering = 0; 





hmod = 
modem.pskmod('M',4, 'PhaseOffset',pi/4,'SymbolOrder','Gray'); 
hdemod = modem.pskdemod (hmod) ; 





while and(or(totalErrors < 5000,NumTests < 400),NumTests <3500) 
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NumTests = NumTests + 1; 
s = fix(rand(1,testSize) *M); % create random signal 
s(l:pilotInterval:testSize) = 3; 


m = modulate (hmod,s)j; 


m = rectpulse(m,fs/Rs) ; 
m = filter(chan,m); 





m = awgn(m,EbNo + 10*1log10(Rs/fs)+ 10*1log10(k)); 





intdump (m, fs/Rs) ; 
r = DD_DF_demod(m, hmod, hdemod, s(1),pilotInterval) ; 


3 
| 
H 


totalErrors = totalErrors + biterr(r,s); 
currentPerr = totalErrors/ (NumTests* (testSize- 

















floor (testSize/pilotInterval) ) *k); 


end 
Pb = currentPerr; 


DD_DF_QPSK_plot.m 





%& Jeff Smith 
%& NPS 
% Creates EbNo plots for QPSK 





ea 





bNo = 10:3:61; 





err = zeros(1,length(EbNo) ); 


for n = 1:length(EbNo) 





err(n) = DD_DF_QPSK(EbNo(n),100) 





end 


semilogy (EbNo(1:length(err)),err) 
hold on; 
semilogy (EbNo,berfading(EbNo, 'psk',4,1),'--'); 




















DD_DF_QAM.m 


function Pb = D_DF_QAM(EDNo, pilotInterval) 








$ D_DF_QAM 
& Jeff Smith 
%& NPS ECE 











% April 23, 2008 


Rs = 1le5; % symbol rate 
fl = 1*Rs; % f1 

fs = f1*1; % sample rate 

fd = 120; 

M = 64; % 


k = log(M)/log(2); 
testSize = 10000; 
currentPerr = 0; 
err = []; 
NumTests = 0; 
totalErrors = 0; 





chan = rayleighchan(1/(fs),fd); 
chan.StoreHistory = 0; 
chan.ResetBeforeFiltering = 1; 





hmod = modem.gqammod('M' 
hdemod = modem.gqamdemod (hmod) ; 


z = 0:M-1; 
y = modulate (hmod, z); 
Pavg = sum(abs(y) .*2)/M; 


[a, b] = max(abs(y)); 





while and(or(totalErrors < 4000,NumTests < 300),NumTests <2500) 





NumTests = NumTests + 1; 
s = fix(rand(1,testSize) *M); 


»M, 'SymbolOrder', 'Gray'); 


% create random signal 





s(l:pilotInterval:testSize) 
m = modulate (hmod,s)j; 


m = rectpulse(m,fs/Rs) ; 
= filter(chan,m); 





,LO*logl10 (Pavg)); 








m = intdump(m, fs/Rs); 


b-1; 


m = awgn(m,EbNo + 10*1l0g10(fs/Rs) +10*10g10 (k) 


r = DD_DF_demod(m, hmod, hdemod, s(1),pilotInterval) ; 


totalErrors = totalErrors 














+ biterr(r,s); 


currentPerr = totalErrors/ (NumTests* (testSize- 





floor (testSize/pilotInterval) ) *k); 
end 
Pb = currentPerr; 


O. DD_DF_QAM_plot.m 


SJeff Smith 
SNPS 
SPlots QAM BER curve 








10:3:61; 





zeros(1,length(EbNo) ); 
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for n = 1:length(EbNo) 


























err(n) = DD_DF_QAM(EbNo (n) ,5) 
end 
semilogy (EbNo (1: length (err) ),err) 
hold on; 
semilogy (EbNo, berfading (EbNo, 'qam',64,1),'--'); 
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